Hadron Formation in Relativistic Nuclear Collisions and the QCD Phase Diagram 
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We analyze hadrochemical freeze-out in central Pb+Pb collisions at CERN SPS and LHC energies. 
We determine the effects of baryon and antibaryon annihilation and/or regeneration occurring during 
the final cascade expansion stage of the collisions by deriving "survival factors" for each hadronic 
species and all energies considered by employing the UrQMD hybrid model. These survival factors 
are shown to resemble the pattern of data deviation from the statistical equilibrium calculations 
with the statistical model. We apply them in the SM data analysis, obtaining new points in the 
T — fiB plane which follow closely the parton-hadron phase boundary predicted by lattice QCD. 
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The phases, and phase transformations of strongly in- 
teracting matter represent one of the key remaining ques- 
tions of the Standard Model. It is the goal of Quantum 
Chromodynamics (QCD) theory to delineate a phase di- 
agram of such matter [l[ . As its most prominent feature, 
recent results of lattice QCD calculations |[2h^] predict 
a phase transformation between confined hadrons and 
deconfined quarks and gluons. A parton-hadron coexis- 
tence line results, in the plane spanned by temperature 
T and baryochemical potential /is, the principal vari- 
ables of a phase diagram derived from the grand canon- 
ical equilibrium thermodynamics of quarks and gluons, 
as considered on the lattice 0]. The coexistence line (or 
phase boundary) originates, at = MeV, with a tem- 
perature T = 165 ± 8 MeV, far into the nonperturbative 
sector of QCD. The nature of the transition is a narrow 
cross-over here 0]. It continues, with a slight downward 
curvature, up to a i^b of about 600 MeV Q perhaps 
featuring a critical point d, Q whereupon the transition 
would become first order. 

Relativistic nucleus- nucleus collisions aim to identify 
such features of the phase diagram [l3| . The large col- 
lisional volume undergoes an evolution of the contained 
QCD matter, starting from conditions far from quark- 
gluon equilibrium during interpenetration of the colli- 
sion partners. After a certain formation time the col- 
lisional fireball will approach quark and gluon chemical 
equilibrium, at least locally, and a hydrodynamic expan- 
sion evolution will set in which proceeds along a trajec- 
tory in the (T, /i^) plane llj. With increasing collision 
energy this trajectory samples across this plane toward 
fiB = Q MeV, the site of the primordial cosmological evo- 
lution, which is closely approached by recent experiments 



various regions of the phase diagram 10]. In this Let- 
ter we shall focus on results concerning the position, in 
{T,ij,b), of the parton-hadron phase boundary line. In 
fact we will show that hadron production data confirm 
the QCD results [1, 0| for this line. 

Hadrons get formed once the expansive evolution 
crosses the phase boundary. Hadronization seems to be 
close to the chemical freczc-out point, that is the point 
where subsequent inelastic collisions between hadrons 
cease and hadronic species abundances get frozen [13, l3|- 
Moreover, the resulting hadronic yield distributions, over 
the various species, can be understood to closely resemble 
grand canonical Gibbs equilibrium ensembles, from AGS 
to LHC energies. This observation of equilibrium 15;] is 
exploited in the Statistical Hadronization Model (SHM) 
[r6-18] which yields a "freeze-out point" in the (T, /is) 
plane, for each collision geometry and energy studied. 
Such points are then smoothly interpolated yielding the 
hadro-chcmical freeze-out curve 19] which is frequently 
shown in the QCD phase diagram. 

The freeze-out curve converges towards the lattice 
QCD phase boundary line, at small /is, thus confirm- 
ing the transition temperature of about 165 MeV [l^ . 
But it appears to fall well below the line toward higher 
values of ^b, a widely discussed feature which has, 
however, not been conclusively understood. 



at RHIC and LHC [I4]. Various physics observables get 



formed at different stages of the evolution. They "freeze 
in" thus surviving the subsequent stages essentially un- 
obliterated, and thus preserve information pertinent to 



We have shown in a previous publication [21j that the 
hadronic freeze-out curve indeed needs revision. All pre- 
vious determinations of its points in the (T, ^b) plane, in 
the framework of SHM, have implicitly assumed that a 
primordial hadro-chemical equilibrium (an intrinsic fea- 
ture of the hadronization process as observed in ele- 
mentary collisions [22| ]) remains frozen-in throughout the 
final expansion phase. In other words, the chemical 
freeze-out point was assumed, in these analyses, to co- 
incide with the point of the latest chemical equilibrium 
for the hadronic species. This assumption turns out. 
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on the one hand, to be reahstic as far as the inelas- 
tic sector during the hadronic phase is concerned, that 
dehvers the bulk mesonic output. On the other hand, 
the baryon-antibaryon annihilation and regeneration pro- 
cesses do not fall away with the onset of expansive cooling 

Hill. 
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Their final effect consists of a considerable 
distortion of the initial, post-hadronization equilibrium 
yield distribution, in the antibaryon and baryon sector 
[25|. We have shown 2lj that such distortions univer- 
sally lead to a split between the latest hadrochemical 
equilibrium point (LHCE) and the chemical freeze-out 
point with a downward shift of the latter [l9| , and to un- 
satisfactory SHM fits. In fact these effects have recently 
also been noticed with the first LHC hadron production 
data of ALICE 2^ 27 1, which resulted in recognition of 
a "non-thermal proton to pion production ratio" (28| . 

We shall demonstrate in this Letter that an appropri- 
ate second order correction of the observed antibaryon 
and baryon yields will lead to a revised LHCE curve 
as derived from SM analysis of data gathered from SPS 



[29| to current LHC [27| energies. We employ correc- 
tion factors derived from analysis of the cascade phase 
effects employing the recent hybrid version of the micro- 
scopic transport model UrQMD (s^, |3l| , that agree, at 
the SPS energy, with results of a former study by Rapp 
and Shuryak 24] that employed a blast wave expansion 
model including detailed balance of baryon-antibaryon 
interaction. The resulting LHCE points will be shown to 
coincide with the recent predictions ^, J\ of the parton- 
hadron boundary line from lattice QCD at finite bary- 
ochemical potential. 

We illustrate in Fig.l the effects of the final 
hadron/resonance expansion phase on the observed yield 
distributions, as derived from the UrQMD hybrid ver- 
sion 3l|. It features a 3-1-1 dimensional hydrodynamic 



expansion during the high density stage, terminated by 
the Cooper-Frye hadron formation mechanism (assum- 
ing local thermodynamical equilibrium) once the energy 
density of flow cells falls below a "critical" value, as- 
sumed here to be Cc = 0.8 GeV/fm^. This criterion mod- 
els the effect of crossing the QCD phase boundary. The 
hadron distribution can be examined at this stage, emit- 
ting into vacuum. Alternatively, the UrQMD cascade 
expansion stage is attached to the Cooper-Frye output, 
as an " afterburner" . The effect of this stage can be quan- 
tified by modification factors, for each hadron multiplic- 
ity, M = N{Hydro + Aft)/ N (Hydro). These factors 
are shown in Fig.l (top panel). It illustrates the results 
obtained for central Pb-|-Pb collisions at SPS energies, 
= 17.3 GeV,8.7 GeV and 7.6 GeV and for 

the present top LHC energy, =2.76 TeV At 

the SPS energies we see the bulk hadron output relatively 
unaffected by the afterburner, including the S, f2 and f2 
yields. Whereas the other antibaryons, p, A and S are 
showing significant suppressions ranging from 50 to 25%. 
At LHC energy the modification factors of baryons and 
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FIG. 1: Top panel: Modification factors from UrQMD be- 
tween particle multiplicities at hadronization and after the 
hadronic cascade afterburner stage for Pb-Pb collisions at 
^/s^ = 2.7 TeV, 17.3 GeV, 8.7 GeV and 7.6 GeV. Bot- 
tom panel; Ratio between central Pb-Pb collisions data at 
^Snn = 17.3 GeV and a statistical model fit excluding an- 
tibaryons [2ll ]. for which the effect of the hadronic stage is 
largest. 



antibaryons in Fig.l become approximately equal as is to 
be expected in view of the particle-antiparticle symmetry 
prevailing here (with fiB close to zero). The suppression 
pattern differs, in detail, from the pattern at SPS en- 
ergy. It appears to be restricted to the p, p, S and 5 
yields whereas the A, D, and their antiparticles exhibit 
influences of a possible dynamical regeneration (see ref 
(23| for discussion). 

Within the above model considerations the annihila- 
tion and/or regeneration effects inflict distortions of the 
initial equilibrium yield distributions imprinted into the 
subsequent cascade phase by the grand canonical Cooper- 
Frye formalism. It is important to demonstrate that a 
quantitatively similar distortion pattern governs the ex- 
perimentally observed hadron multiplicity data. To this 
end we have performed a SHM analysis of the NA49 
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hadron yield data [29[ for central Pb+Pb collisions at 
17.3 GeV excluding the most affected antibaryon species 
p, A and S from the fit procedure (see ref [21| for de- 
tail). The result is shown in the bottom panel of Fig.l. 
It shows the ratios of data relative to SM predictions, 
for all species. We note that the bulk hadron data are 
well reproduced whereas the yields of p, A and 2 ex- 
hibit a strong suppression relative to the SM equilibrium 
distribution of multiplicities. The pattern quite closely 
resembles the UrQMD prediction in the upper panel of 
Fig.l. 

These observations lead to the idea to employ the 
UrQMD "survival factors" on face value: as correction 
factors employed in the SHM analysis that aims at con- 
structing a LHCE curve. Such an analysis is shown in 
Fig. 2, applied to recent LHC ALICE data [27] for the 
20% most central Pb-HPb coUisions at 2.76 TcV. The 
top panel gives the standard SHM fit which is unsatis- 
factory, the bulk pion and kaon yields being missed at the 
cost of accounting for the baryon sector. Similar results 
have been obtained in ref. [s^l • The bottom panel shows 
the analysis with correction factors (the survival factors 
from UrQMD in Fig.l) applied to the SHM fit procedure. 
It yields a (T, hb) point at (166 MeV, 2 MeV), with im- 
proved x^- 

We repeated this analysis with the NA49 data for 
5% centrality selected Pb-|-Pb collisions [i^ at y^Snn = 
7.6,8.7 and 17.3 GeV. All obtained SM parameters are 
gathered in Tab.l. We presently have not used the RHIC 
data at 200 GeV because of difficulties in combining 
measurements from different experiments. More data 
are forthcoming from the RHIC beam energy scan pro- 
gram [31] which can be used to systematically extend the 
present SHM analysis in the future. 

Fig. 3 shows our principal result. The four obtained 
LHCE points are inserted into a phase diagram obtained 
recently by QCD lattice calculations at finite baryochem- 
ical potential, by Enroedi et al.Q. Similar results can be 
found in . The authors distinguish two different deter- 
minations of the position of the critical curve, based on 
the chiral condensate (^\1/) and on the strange quark 
susceptibility Xs/T^, respectively. The resulting two 
close curves Tq (pb) cover baryochemical potentials up 
to 600 MeV. The LHCE points follow the latter theoreti- 
cal choice for the parton-hadron transition line. They are 
listed in Tab.l. One generally observes an upward shift 
of T from chemical freeze-out to LHCE upon application 
of the UrQMD correction factors. Except at 8.7 GeV it 
is accompanied by an improved x^. 

A brief reflection is in order here concerning our em- 
ploy of the UrQMD hybrid transport model [sij. Its 
account for the final hadron/resonance cascade evolution 
does not include the reverse of the annihilation processes, 
such as B + pio 5 pions, that could modify the survival 
factors 2J, |3J|. This is a general feature of all exist- 
ing microscopic transport models. To account for such 
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FIG. 2: Statistical model fits to preliminary ALICE data [23| 
for 20% central Pb-Pb collisions at -^Snn = 2.7 TeV (a) and 
to the same data but with modification factors from UrQMD 
applied in the statistical model fits (b). 
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FIG. 3: Phase diagram of nuclear matter in the (T, fis) plane 
with predictions from Lattice QCD calculations (solid line: 
based on strange quark susceptibility xs, dashed line: based 
on chiral condensate {^'i')) 5, 6] and from statistical model 
fits at SPS and RHIC energies (dotted line) 18]. The results 
of the modified statistical model fits of this letter are shown 
as closed circles. 



reverse reaction channels and establish, in principle, ef- 
fects of detailed balance, one has to employ analytic mod- 
els of fireball expansion \24l |3^. These, in turn, imply 
thermodynamically homogeneous collision volumes that 
miss the effects of local and surface density fluctuations. 
They thus overestimate both the annihilation and, in par- 
ticular, the regeneration rates which scale with the fifth 
power of the density. Most remarkably, the main predic- 
tion 13, 13] is, again, a net loss of about 50% in the p 
yield, similar to Fig.l. 

In summary we have demonstrated that the semi- 
empirical freeze-out curve from SHM analysis of the 
hadronic multiplicity data at various incident energies 
does not coincide with the latest point at which hadrons 
are in chemical equilibrium following the parton-hadron 
conversion. To obtain the latter, one has to correct 
for abundance changes originating in the course of the 
hadron/resonance expansion that ends the dynamical 
evolution of A-l-A collisions. We have quantified those 
changes by survival factors obtained from the UrQMD 
transport model. Significant modifications are restricted 
to the baryon-antibaryon sector, and their detailed pat- 
tern depends on the incident energy. We have chosen 
to take account of the cascade phase effects by calcu- 
lating the UrQMD survival factors and employing them 
to correct the SM predictions. At the four energies 
considered here the resulting freeze-out points coincide 
with the lattice QCD phase boundary line, thus resolv- 
ing a long standing problem. We cover the interval 
< mub < 430 MeV. It will be very interesting to 
continue this analysis further upward in /is (and thus 
downward in incident energy). 
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Std. fit 
Mod. fit 


156 ±5 
166 ±3 


1 ± 12 
2± 6 


1.09 ±0.07 
0.98 ±0.04 


26.5/9 
11.5/9 


Pb-Pb 5% central = 17.3 GeV 


Std. fit 
Mod. fit 


151 ±4 
163 ±4 


266 ±9 
250 ±9 


0.91 ±0.05 
0.83 ±0.04 


26.9/11 
20.4/11 


Pb-Pb 5% central ^Ts^^ = 8.7 GoV 


Std. fit 
Mod. fit 


148 ±4 
161 ±6 


385 ± 11 
376 ± 15 


0.78 ±0.06 
0.72 ±0.06 


17.9/9 
25.9/9 


Pb-Pb 5% central ^i^^^ = 7.6 GeV 


Std. fit 
Mod. fit 


140 ± 1 
156 ±5 


437 ±5 
426 ±4 


0.91 ±0.01 
0.81 ±0.00 


22.4/7 
14.7/7 



TABLE 1: 
SPS data. 
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